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Abstract 

We consider a Langevin equation with variable drift and diffusion coefficients separable in time 

O ' and space and its corresponding Fokker-Planck equation in the Stratonovich approach. From 

^ ! 

^ ' this Fokker-Planck equation we obtain a class of exact solutions with the same spatial drift and 

C^ • diffusion coefficients. Furthermore, we analyze some details of this system by using the spatial 
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diffusion coefficient D(x) = \^\x\ 2. 
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In these last decades, anomalous diffusion properties have been extensively investigated 
by several approaches in order to model different kinds of probability distributions such 
as long-range spatial or temporal correlations [l|, |2|. These approaches have been used 
to describe numerous systems in several contexts such as physics, hydrology, chemistry and 
biology. The diffusion process is classified according to the mean square displacement (MSD) 

(x^(t)>~r. (1) 

In the case of normal diffusion, the MSD grows linearly with time (a = 1). For < a < 1 
the process is called subdiffusive, and for a > 1 the process is called superdiffusive. The well- 
established property of the normal diffusion described by the Gaussian distribution can be 
obtained by the usual Fokker-Planck equation with a constant diffusion coefficient (without 
the drift term) |3|, \^ or by an integro-differential diffusion equation with the exponential 
function for the waiting time probability distribution |5|. Anomalous diffusion regimes can 
also be obtained by the usual Fokker-Planck equation, however, they arise from variable 
diffusion coefficient which depends on time and/or space. On the other hand, in the view 
of Langevin approach it is associated with a multiplicative noise term. In other approaches 
such as the generalized Fokker-Planck equation (nonlinear) and fractional equations, they 
can describe anomalous diffusion regimes with a constant diffusion coefficient. 

The Langevin equation is a very important tool for describing systems out of equilibrium 
[sl, |j]. Moreover, this equation has been extensively investigated; many properties and 
analytical solutions of it have also been revealed. In this work, we present solutions of a 
class of the Langevin equation with the deterministic drift and multiplicative noise terms 
in time and space. To do so, we obtain the corresponding Fokker-Planck equation in the 
Stratonovich definition, and then we obtain its solutions for the probability distribution 
function (PDF). 

Langevin equation and its corresponding Fokker-Planck equation. We consider the follow- 
ing Langevin equation in one-dimensional space with a multiplicative noise term: 

i = hi^,t) + g{^,t)Tit), (2) 

where ^ is a stochastic variable and r(t) is the Langevin force. We assume that the averages 
(r(t)) = and (r(t)r(t)) = 25{t-t) [3]. h{^,t) is the deterministic drift. Physically, the 
additive noise {g{C,,t) constant) may represent the heat bath acting on the particle of the 



system, and the multiplicative noise term, for variable g{^,t), may represent a fluctuating 
barrier. For g = \AD and h{^,t) = 0, Eq. ([2]) describes the Wiener process and the 
corresponding probability distribution is described by a Gaussian function. In the case of 
g{^,t), some specific functions have been employed to study, for instance, turbulent flows 

g{x,t) ~ |x|°t'') la-lSj. By applying the Stratonovich approach in a one-dimensional space 

3(], we obtain the following dynamic equation for the PDF: 



dW{x,t) 
di 



-— [D,{x, t)W{x, t)] + — [D^ix, t)W{x, t)] 



where Di{x,t) and D2{x,t) are the drift and diffusion coefficients given by 

Di{x,t) = h{x,t) + ^^^] g{x,t) 



dx 



and 



(3) 



(4) 



(5) 



D2{x,t) = g^{x,t). 

We note that Eq. (^ has a spurious drift due to the Stratonovich definition. Moreover, 
Eq. ([3]) can be written as 



dW{x,t) 
di 



-^[hix,t)Wix,t)] + ^ 



g{x,t) 



dg{x,t)W{x,t) 
dx 



(6) 



i^'or the case of h{x,t) = and g{x,t) = T(t)D{x) the system has been considered in Ref. 
9|; the solution for W{x,t) is given by 



W{x,t) = B{t) 



exp 



x{xY 
' 4t{t) 



D(x)vTity 



where 



dt 
It " 
dx 
dx 



T\t) , 

1 

" D{x) 



(7) 

(8) 
(9) 



and B (t) is a normalization factor. Eq. ([71) can describe interesting properties such as 
non-Gaussian distribution, combination of behaviors like Gaussian (for small distance) and 
exponential (for large distance), and combination of behaviors like Gaussian (for small dis- 
tance) and power law decay for long distance. Further, it can describe many bimodal 

/ — -- 
distributions for different forms of g{x,t). For instance, let us consider D{x) = \/ D \x\ 2, 



then the probabihty distribution and MSD are given by 

._ r N^ 

X2- 

VAnDt 
and 



wixrt) = \x\^-^-^^^^^rii (10) 



[D' (2 + 6)'] '+" r (4^^ t^ (t) 
[At)) = L^ ^ ^^'''^''^ -, ^ > -2 . (11) 



Moreover, the PDF can also be obtained for 6 = —2; in this case the PDF gives a log-normal 

/ — — - 
distribution. One can see that the multiplicative noise term in space D{x) = \/ D \x\ ^ 

produces non-Gaussian shapes for the PDF (TTOl) : it presents a Gaussian shape only for 

^ = 0. It can also reproduce the asymptotic behavior of the random- walk model and time 

fractional dynamic equation for t = t''(2+6»)/2 jgj^ where < /3 < 1. Now we want to show 

two interesting processes which can be obtained from Eqs. (TTOj) and (TTTj) . To do so, we take 

9 > —2. The first one we consider a simple expression for T(t) given by 



T{t) = ^, (12) 



for t ^ 1 . From Eq. ([8]) yields 



t{t) = q\nt. (13) 

Eqs. (fTTl) and (fT3l) describe the ultraslow diffusion processes. This kind of diffusion has 



been found, for instance, in aperiodic environments 10|. 
The second one we consider the following T{t): 






Vat" \/E"=oCjAj 
nt) = yl ,-v- , (14) 

where Cj, \j and a are constants. Using the function f lT4|) one can obtain anomalous diffusion 
processes with logarithmic oscillations. We note that the time behavior with logarithmic 
oscillation is ubiquitous: examples have been observed, for instance, in epidemic spreading 



in fractal media [ll|, financial stock market [12] and diffusion-limited aggregates [13]. In 
Fig. 1 we show the function T{t) flT4l) for Aj = a\ Ci = {a/by, a = 1/15 and b = 0.3; 
for these values the curves present logarithmic oscillations with different values of n and a. 
From Eq. ([8]) we obtain 

t{t)= n I ,^,. ■ (15) 



Moreover, the PDF (fTOj) presents unimodal states for —2 < 6 < and bimodal states 
for 6 > (see in Fig. 2) with pronounced cusps. The numerical results show that the PDF 
changes practically nothing for n = 2 and n = 6. In Fig. 3 we show the MSD f llip in function 
of time t] it presents anomalous diffusion processes with logarithmic oscillations. We see that 
the main trends in the MSD have power-law behaviors which indicate subdiffusive regimes. 

We now consider that the deterministic drift h{x,t) and multiplicative noise term g{x,t) 
are separable in time and space, and they are given by 



and 

Then, Eq. (|3]) reduces to 
dW{x,t) 



h{x, t) = Ti (t) D (x) 



g{x,t)=T{t)D{x). 



d 



d 



dt 



Ti{t)^[D{x)W{x,t)]+T\t)^ 



Dix 



dD{x)W{x,t) 
dx 



(16) 



(17) 



(18) 



We note that the coefficients h{x,t) and g{x,t) given by h{x,t) = g{x,t) = D(x) have 
been used for studying Brownian pumping in nonequilibrium transport processes IJ]. By 
suitable transformations of variables we can show that Eq. ( [T8l) can be reduced to the 
const ant- diffusion equation without the drift coefficient term. To do so, we take the following 
transformations: 



p(x,t) = D{x)W{x,t) 
dt 



T^t) 



and 



X 



dx 
Dix) 



- / dtTi (t) + A, 



(19) 
(20) 

(21) 



where A is a constant, then Eq. (ITSl) reduces to 



dp{t*,x*) d'^p{t*,x*) 



(22) 



dt* dx*"^ 

Eqs. (1201) and ( 12T]) give the time and space scaling factors which connect Eq. (ITSl) to the 

ordinary diffusion equation ( 122|) . Eq. ( l22l) can be solved and the solution with a natural 

boundary condition is given by 



p{t*,x*) =C- 



exp 



X 



Vt* 



(23) 



where C is a normalization factor. Eqs. (fT9|) and ( 123|) show that the time-dependent 
coefficients T[t) and Ti(t) do not change the Gaussian form, however the coefficient D{x) 
can produce different forms for the distribution W{x,t) 9|]. We note that for D(x) = y/D, 
T{t) = 1 and Ti{t) = the Wiener process is recovered. 

_ 9 

In order to investigate some details of the solution ( 1231) we take D{x) = \j D \x\ ^ . From 
Eqs. ([T9D and dnD, with A = 0, yields 

-a; 2 



W{x,t) 



W{x,t) 



'A-KDt* it) 



X2 



^AixDt* (t) 



exp 



D{2 + e)2t* (t) 



X < , 



exp 



X 2 — 



v^(2+6>) 



/7(t) 



21 



D(2 + e)2t* (t) 



X > 



(24) 



where H{t) = JdtTi (t). Eq. fl2^ shows that the drift term produces an asymmetric PDF 
with respect to the coordinate x. For Ti (t) = the PDF fl24p reduces to the solution fllUp 
without the presence of the drift term, and the symmetric PDF is recovered. In this case, 
the drift term Ti (t) gives the duration of this asymmetry. In Fig. 4 we show the asymmetric 
PDF Eq. (^^ for t = 0.2. The asymmetry of the PDF with 6 = —0.1 is more pronounced 
than the PDF with 6 = -0.5. From Eq. ([24]) we obtain 






m 



6+e 

2(2+6»; 



TT 



L ^ ^ ^ '^ ^ ^ V2(2 + ^) 2' t*(t) 



(25) 



where i-Fi(a, 6, z) is the Kummer confluent hypergeometric function 15|. For Jiif) = 0, 
without the drift term, we recover the result Eq. flTTl) . Moreover, Eqs. fl24|) and fl25|) 
also present an interesting result; for H'^{t)/t* (t) proportional to a constant they give the 
similar results of Eqs. (fTOj) and ( fTTj) . without the drift term. In this case, the drift term 
only contributes an additional constant to the overall behavior of the system. 

We should note that the solutions fITIJl) and (1251) can adequately work for D{x) positive. 
For D{x) containing negative values we should modify and take p {x, t) = —D{x)W{x, t) for 
D{x) negative. For example, let us consider D{x) = x. Then we take 



p (x, t) = -xW{x, t), X <0 



(26) 



and 



p (x, t) = xW{x, t), X > 



(27) 



From Eq. (12T!) we obtain 

X* = ln\x\-H{t)-\n\xo\ (28) 

and 

{\n\x\-H{t)-ln\xo\f 



exp 



44* (t) 



W{x,t) = ^ ^. (29) 



4A/7rt* (t) \x\ 
This is the log-normal distribution. The distribution fl2^ is the same as the one given in 



Ref. [16|| for t* (t) = t, which has been obtained from the method of characteristic. It is 
worth mentioning that the coefficients h{x, t) and g{x, t) given by h{x, t) ~ a; and g{x, t) r^ x 
might be used to investigate the barrier crossing problem in heavy- ion fusion reaction jlTj ]. 
and they can also be used as a limiting case of a Langevin equation for describing the tumor 
cell growth system |18 |. 

Conclusion. When a multiplicative noise term is introduced into the simple Langevin 
equation (|5]), even separable in time and space, the system can exhibit complex behaviors 
and a rich variety of processes. We have analytically presented a class of these processes. 
We hope that they can be used to mimic a wide class of natural systems. 
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FIGURE CAPTIONS 



FIG. 1 - Plots of the function T(t), Eq. flT^ . The dashed hnes correspond to a = 0.5, 
whereas the sohd hnes correspond to a = 1. 

FIG. 2 - Plots of the PDF ([10]) for A^ = a\ d = {a/h)\ a = 1/15, b = 0.3, D = 1,9 = 0.5 
and a = 1. The solid lines correspond to n = 2, whereas the dotted lines correspond to 
n = 6. 

FIG. 3 - Plots of the MSD (^ for A, = a\ c, = {a/h)\ a = 1/15, b = 0.3, D = 1 and 
9 = 0.5. The solid line with n = 4 corresponds to a = 0.5, whereas the solid line with n = 6 
corresponds to a = 1. The dashed lines correspond to the power-law functions. 

FIG. 4 - Plots of the PDF (^ for D = 1, t* (t) = t and H{t) = t. The dotted line 
corresponds to ^ = —0.5, whereas the solid line corresponds to ^ = —0.1. 
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